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^ '. Abstract 
O 

0^ \ A Monte Carlo study of the azimuthal angular distribution around the virtual 

boson-proton beam axis for dijet events in DIS at HERA is presented. In the 
presence of typical acceptance cuts on the jets in the laboratory frame, the 
Oh. azimuthal distribution is dominated by kinematic effects rather than the typi- 

(~| \ cal cos (j) and cos 2^ dependence predicted by the QCD matrix elements. This 

implies that the (j) dependent part of the QCD matrix elements contributes 
■ even to the dijet production cross section. Neglecting this (j) dependence leads 

j_j ' to an error of about 5-8% in the production cross section for typical accep- 



tance cuts in the laboratory frame. We also present first NLO results on the 
(^-decorrelation of the jets through NLO effects. 



Jet production in deep inelastic scattering (DIS) is an important laboratory for testing 
our understanding of perturbative QCD. Good event statistics allow for precise measure- 
ments and thus for a variety of tests of our understanding of QCD dynamics. For quantitative 
studies one clearly needs to compare data with calculations which include next-to-leading 
order (NLO) QCD corrections. The importance of keeping the full helicity structure in such 
a calculation is demonstrated in this letter for the example of dijet production in DIS. 

Dijet production in neutral current (NC) DIS lepton proton scattering proceeds via the 
exchange of an intermediate vector boson 7*/^. We denote the ^* jZ momentum by g, 
its absolute square by Q^, and use the standard scaling variables x = Q'^/{2P-q) and 
y = P ■ q/ P ■ I. In Born approximation, the subprocesses 

e(/) + g(po) ^ e{V) + q{pi) + g{p2) (1) 
eil) + g{Vo) ^ e(/') + + q^p^) (2) 

contribute to the NC two-jet cross section. NLO effects will be discussed later. In the absence 
of jet cuts in the laboratory frame, the full leading order (LO) QCD matrix elements predict 
a typical azimuthal angular distribution of the jets of the form |jl]J^ 

= A + 5cos0 + Ccos20 (3) 

(icos0 

where denotes the azimuthal angle of the jets around the virtual boson direction. This 
angular distribution is determined by the gauge boson polarization: the coefficients A,B,C 
in Eq. (|^) are linearly related to the nine polarization density elements hmm' of the exchanged 
gauge boson: 

hmm' = el{m)H^''e,{m') (m, m' = +, 0, -) (4) 

where 

e^(±) = i=(0; ±1, -I, 0) e^(0) = (0; 0, 0, 1) (5) 

are the polarization vectors of the exchanged boson in the hadronic [i.e. boson-proton) 
center of mass frame (HCM) and H^^ denotes the hadronic tensor for the contributing QCD 
processes in Eqs. (|1],0). The coefficient A in Eq. (|^) can be expressed in terms of the diagonal 

density matrix elements /loo, ^++5 h , the coefficient B in terms of transverse- longitudinal 

interference matrix elements /i+o, ^o+; ^-o, ^0- and the coefficient C in terms of transverse 

interference density matrix elements h^ , @- Technically, the coefficient A can be 

calculated from the hadronic tensor in the one photon exchange case by the two covariant 
projections g^y and Po^Pov on the hadronic tensor, whereas projections defined with final 
state hadronic momenta are necessary to calculate B and C [@. 

The NLO eP — > n jets event generator MEP JET 0] is used for the subsequent numerical 
studies in this paper. The corresponding cross sections and distributions are calculated by 
evaluating the corresponding one-loop and tree-level helicity amplitudes and therefore, the 
MEP JET program allows for the calculation of all possible jet-jet and jet-lepton correlations 
in NLO, including the (j) dependence in Eq. (^. 
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FIG. 1. a) Normalized (p distribution in LO around the virtual boson direction of jets in the 
cone scheme defined in the HCM with 

Pt^^(j) > 5 GeV. The integrated dijet cross section is 
1465pb; b) same as a) but but with an additional cut of py'^(j) > 5 (sohd) {pT^{j) > 4 (dotted)) 
GeV on the jets. The (p integrated dijet cross section is 945 pb (1146 pb); c) same as a) but with 
an additional |??^'^'^(j)| < 2 cut on the jets. The (j) integrated dijet cross section is 1087 pb. All 
results are given in LO with 7* and Z exchange. 



The following numerical studies for e'^P scattering with 7/Z exchange are based on 
MRSRl 1^ parton distribution functions with the two loop formula for the strong coupling 
constant and Uf = 5. The renormalization and the factorization scales are set to fiR = 
ftp = (^rO)); the average of the jets, which appears to be the "natural" scale for 
jet production in DIS 0. Here, (fcy(j))^ is defined by 2Ej{l — cosOjp) in the Breit frame, 
where the subscripts j and P denote the jet and proton, respectively. The lepton and hadron 
beam energies are 27.5 and 820 GeV, respectively and events are selected in the following 
kinematical region: 40 GeV^ < < 2500 GeV^ 0.04 < y < I, and E{e') > 10 GeV, where 
E{e) is the energy of the scattered lepton. Unless stated otherwise, no cuts on the pseudo- 
rapidity T]^^^ = — lntan(^^/2) or the transverse momentum on the jets in the lab frame are 
imposed. 

Without an (experimental) separation of a quark, anti-quark or gluon jet, the cos(/> term 
in Eq. (|^) is not observable and thus has to be averaged out. The resulting LO cos 20 
dependence of dijet events is shown in Fig. |T|a for a cone scheme defined in the HCM with 
radius R = 1 and Pt^^{]) > 5 GeV. We find very similar results for jets in the kp algorithm 
implemented in the Breit frame with = 40 GeV^ and ?/cut = 1- Note that the HCM 
and the Breit frame are related by a boost along the boson-proton direction which does 
not change the dependence in Eq. (|^). One observes a sizable cos 20 dependence (after 
averaging over cos0) in Fig. |T]a as expected from Eq. (^). The size of the dependence 
in this normalized distribution is rather insensitive on the Px^^ (or Ep in the kp scheme) 
requirements on the jets. Averaging (integrating) over implies that the coefficients B and 
C in Eq. (|]) do not contribute to the dijet production cross section in this case. 

Fig. |l]b shows the distribution of the jets defined in the cone scheme and cuts as in 
Fig. la, but with an additional p^^ij) > 5 (sohd) {pt^U) > 4 (dotted)) GeV cut on the 
jets in the laboratory frame. Only 65% (78%) of the events pass the additional Pr^ij) 
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TABLE I. Effects of the jet azimuthal angular distribution on the dijet production cross section 
in DIS as a function of py'^(j) cuts on the jets, cr^^^'^* (^^japprox-j i^^ggj ^^^^ exact (approximate) 
matrix elements including (without) the coefficients B and C in Eq. (^). Jets are defined in the 
HCM with radius R = 1 and and Pr'^^ij) > 5 GeV. Additional parameters are explained in the 
text. Similar results are found for jets defined in the fc^ scheme. 



lab frame cut 



exact 



[2-jet] 



approx 



[2-jet] 



Ad 



Jijapprox _ ^cxactj^^c 



p'fU ) > 4 
Pt^(j) > 8 

p'S^'U) > 10 



1465 pb 
1390 pb 
1145 pb 
689 pb 
348 pb 
187 pb 



1465 pb 
1405 pb 
1195 pb 
733 pb 
370 pb 
198 pb 



% 

1.1 

4.4 

6.4 

6.3 

5.9^ 



TABLE II. Same as table | as a function of |V'^(i)l cuts on the jets. 



lab frame cut 


^exactj2.jet] 


^approx [2-jet] 


— J^approx _ ^exactjy'^exact 


IV'^^J)! < 5 


1465 pb 


1465 pb 


% 


< 3.5 


1437 pb 


1439 pb 


0.15% 


< 2.5 


1282 pb 


1288 pb 


0.5% 


\v'''\j)\ < 1-5 


801 pb 


820 pb 


2.4% 


Ir^i'^^j)! <0.5 


131 pb 


141 pb 


7.6% 



cuts. These cuts have a dramatic effect on the shape of the distribution. The shape 
of the distribution is now governed by the kinematics of the surviving events. The pJ;^^(j) 
cut introduces a strong dependence. The "kinematical" distribution in Fig. |I]b is very 
different from the "dynamical" cos 20 distribution in Fig. [l|a. The only remaining vestiges 
of the gauge boson polarization effects in the distribution are the dips at = 90° and 270° 
in the dotted curve in Fig. |I]b. Imposing cuts on the pseudo-rapidities |?7^'^''(j)| of the jets in 
the laboratory frame has also a large effect on the shape of the distribution as shown in 
Fig. Ic. 

A consequence of this kinematical dependence is, that the dependent coefficients B 
and C in Eq. (^ now also contribute to the dijet production cross section after integration 
over 0. We have analyzed this effect by comparing the dijet cross section using the exact 
matrix elements (including the dependent coefficients in Eq. (^) with the resulting dijet 
cross section where the coefficients B and C in the matrix elements in Eq. (|]) are neglected. 
The results are shown in table | for various Px^ij) cuts on the jets. We find that the dijet 
cross section based on the approximate matrix elements is larger compared to the correct 
result. The error depends on the p^^ij) cuts and can reach about 6.5 %. The error from 
the approximate matrix elements in the presence of pseudo-rapidity cuts on the jets in the 
lab frame is shown in Table y. A combination of Px^{j) and |?7^'*'^(j)| cuts can lead to even 
larger effects. 

Let us now discuss some NLO effects on the dijet azimuthal angular distribution in e^P 
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scattering with 7/Z exchange, where the full helicity structure is kept. The NLO 0{a1) 
2-jet cross section receives contributions from the one-loop corrections to the Born processes 
in Eqs. (|l|) and (0) and from the integration over the unresolved region (defined by a given 
jet algorithm) of the 3-parton final state tree level matrix elements 

e(/) + q{po) e{l') + q{pi) + g{p2) + g{pz) (6) 
e(0 + 9(Po) e(/') + q{pi) + ^(^2) + g^a) (7) 
e(0 + g{po) e(/') + q{p{) + 9(^2) + ^(Ps) (8) 

and the corresponding antiquark processes with q ^ q. 

The NLO corrections introduce a new (numerically small) additional sin0 and sin 20 
dependence in the NC cross section in Eq. (|^) through the imaginary parts of the one-loop 
contributions [0]. 

In addition, the tree level contributions to dijet production in Eqs. (^U) imply that the 
two jets are no longer necessarily back-to-back in (p like in LO, i.e. one expects A0 = 
|0(jet 1) — 0(jet 2)1 7^ 180°. A deviation from A(f) = 180° can arise for example in the cone 
scheme if one of the three final state partons is well separated in AR from the other two 
partons but does not pass the acceptance cuts (like the pf*"^(j) > 5 GeV cut), whereas 
the remaining two partons passes all jet requirements. Thus the event will be accepted as a 
two-jet event where the jets are, however, no longer balanced in 0. 
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FIG. 2. a) A(/) distribution in NLO for jets defined in the kx scheme with = 40 GeV^. 
Similar results are found for jets in the cone scheme defined in the HCM. b) same as a) for the 
JADE scheme and i/cut = 0.02. 

The effect is shown in Fig. |] for the kx scheme defined in the Breit frame with = 40 
GeV^ (Fig. |a) and the JADE scheme |g] defined in the lab frame with y^ut = 0.02 (Fig. ||b). 
In the latter case, jets are required to have a minimum transverse momentum of 2 GeV in 
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the HCM and the lab frame. The corresponding NLO (LO) cross sections are 1350 pb (1240 
pb) in the scheme and 1570 pb (970 pb) in the JADE scheme. 

Fig. 1^ illustrates that the decorrelation effect through the NLO corrections depends 
strongly on the chosen jet algorithm. The decorrelation is larger in the /ct scheme (or a 
cone scheme) than in the JADE scheme. Note that the central bin around 180° in Fig. 0a 
has a negative weight. This shows that the fixed NLO predictions are not infrared save 
for A0 close to 180°. This effect is caused by the negative contributions from the virtual 
corrections, which contribute only to this bin due to the Born kinematics. For the JADE 
scheme the negative contributions in the central bin are already overcompensated by the 
positive tree-level contributions in Eqs. (§[-§). Since the decorrelation effect is larger in 
Fig. one has either to choose wider bins in A0 for the kx scheme to arrive at a positive 
result in the central bin or alternatively one would have to use resummation techniques to 
obtain a reliable perturbation expansion close to A0=18O°. 

Finally, the small asymmetry in the A0 decorrelation in Fig. |^ is caused by our fixed 
ordering of jet 1 and jet 2. The distributions would be perfectly symmetric without a 
separation of a quark, anti-quark or gluon jet. 

To summarize: The azimuthal angular distribution of dijet events around the boson proton 
direction has been discussed. In the presence of cuts on the jets in the lab frame, the angular 
distribution is dominated by kinematic effects and the residual dynamical effects from the 
gauge boson polarization are small. The kinematic effects imply that the non-diagonal 
helicity density matrix elements of the exchanged boson, which are fully taken into account 
in the NLO event generator MEPJET, contribute also to the dijet production cross section. 
First studies of the dijet angular decorrelation through NLO corrections show a fairly strong 
dependence on the chosen jet algorithm. 
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